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We consider the phonon contribution to the bulk viscosities C,i , (2 and (3 of a cold relativistic 
superfluid. We assume the low temperature T regime and that the transport properties of the 
system are dominated by the phonons. We use kinetic theory in the relaxation time approximation 
and the low energy effective field theory of the corresponding system. The parametric dependence 
of the bulk viscosity coefficients is fixed once the equation of state is specified, and the phonon 
dispersion law to cubic order in momentum is known. We first present a general discussion, valid 
for any superfiuid, then we focus on the color-fiavor locked superfluid because all the parameters 
needed in the analysis can be computed in the high density limit of QCD, and also because of the 
possible astrophysical applications. For the three independent bulk viscosity coefficients we find 
that they scale with the temperature as ~ 1/T, and that in the conformal limit only the third 
coefficient (^3 is non-zero. 

. PACS numbers: 47.37.+q,97.60. Jd,21.65.Qr 
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I. INTRODUCTION 



' Superfluidity was first discovered in "^He when this was cooled at temperatures below 2.17 K [H, [Ij. Later on it was 
found that various quantum liquids, both bosonic and fermionic, exhibit the same property. The phenomenon is due 
^ • to the appearance of a quantum condensate, which spontaneously breaks a global U{1) symmetry. 
1^ I Relativistic superfluids might be realized in the interior of neutron stars where the temperature is low and the energy 
i-S^ . scale of the particles is sufficiently high. In particular, in the inner crust of neutron stars the attractive interaction 
between neutrons can lead to the formation of a BCS condensate, the system then becomes superfluid. Moreover, 
. if deconfined quark matter is present in the core of neutron stars it will very likely be in a color superconducting 
^ ' phase 0. Quantum Chromodynamics (QCD) predicts that at asymptotically high densities quark matter is in the 
^sD ] color-flavor locked phase (CFL) Q. In this phase up, down and strange quarks of all three colors pair forming a 
QQ ' difermion condensate that is antisymmetric in color and flavor indices. The order parameter breaks the baryonic 
number U{1)b symmetry spontaneously, and therefore CFL quark matter is a superfluid as well. 

If superfluidity occurs in the interior of compact stars, it should be possible to find signatures of its presence through 
On ' a variety of astrophysical phenomena. For example, the most natural explanation for the sudden spin- up of pulsars [5| , 
the so-called glitches, relies on the existence of a superfluid component in the interior of the star, rotating much faster 
than the outer solid crust. After the unpinning of the superfluid vortices, there is a transfer of angular momentum 
from the interior of the star to the outer crust, giving rise to the the pulsar glitch. 

Another possibility to detect or discard the presence of relativistic superfluid phases consists in studying the 
evolution of the r-mode oscillations of compact stars @. R- modes are non-radial oscillations of the star with the 
Coriolis force acting as the restoring force. They provide a severe limitation on the rotation frequency of the star 
d ' through coupling to gravitational radiation (GR) emission. When dissipativc phenomena damp these r-modes the 
star can rotate without losing angular momentum to GR. If dissipativc phenomena are not strong enough, these 
oscillations will grow exponentially and the star will keep slowing down until some dissipation mechanism is able to 
damp the r-mode oscillations. Therefore, the study of r-modes is useful in constraining the stellar structure and can 
be used to rule out some specific matter phases. For such studies it is necessary to consider in detail all the dissipativc 
processes and to compute the corresponding transport coefficients. 

On the theoretical ground, there is renewed interest in the topic of this manuscript, as it has been recently proposed 
an holographic model to describe relativistic superfiuidity 0, • With the holographic techniques one intends to model 
a strongly interacting gauge theory with a weakly interacting gravity dual, having in this way a prescription to study 
many properties of the gauge system, such as its transport coefficients. 

There are different formulations of the hydrodynamical equations governing a relativistic superfluid when there is 
no dissipation jol-fl^. They were derived as relativistic generalizations of Landau's two-fluid model of non-relativistic 
superfluid dynamics [U, In the non-dissipative limit it is possible to show that all these approaches are equivalent. 
The dissipativc terms which enter into the relativistic hydrodynamical equations of one the above approaches have 
been derived in Ref. [T5| . As it occurs in the non-relativistic case, for a relativistic superfluid one can define a thermal 
conductivity, k, a shear viscosity, ry, and four bulk viscosity coefficients, Cij C2j Csj C4- While C2 has the same meaning 
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as the bulk viscosity of a normal fluid, Ci, Ca ^^'^ C4 refer to dissipative process that lead to entropy production only in 
the presence of a space-time dependent relative motion between the superfluid and the normal fluid components [l|, [2| . 

This paper is devoted to determine the phonon contribution to the bulk viscosity coefficients of a cold relativistic 
superfluid. The superfluid phonon is the Goldstone mode associated with the spontaneous breaking of a U{1) symme- 
try. Being a massless mode, at low temperatures it gives the leading contribution to both the thermodynamics and 
transport phenomena. The contribution of other quasiparticle degrees of freedom, which typically have an associated 
non-zero energy gap E, is suppressed when T <^ E and we will assume that we are in such a temperature regime. 
This work is somehow a follow up of Ref. [l7j . where similar computations were done for a non- relativistic superfluid, 
focusing on ultracold Fermi gases close to the unitary limit. 

We perform our computation in the framework of kinetic theory, with the use of the relaxation time approximation 
(RTA). We present a general discussion regarding the phonon contribution to the transport coefficients for a generic 
relativistic superfluid. The basic ingredient in our derivation is the low energy effective field theory associated to 
the superfluid phonon, which is essentially dictated by symmetry considerations [H, [l^ . The effective Lagrangian 
is presented as the typical expansion in derivatives of the Goldstone field, and with it one can derive the phonon 
dispersion law and all the scattering rates that are needed for the computation of the transport coefficients to a given 
accuracy. It turns out that the scaling laws for the transport coefficients depend on the coefficients of the effective 
field theory, which have to be matched with the microscopic theory. 

After the presentation of the discussion of the bulk viscosities for a generic relativistic superfluid, in the last part 
of the present paper we focus on the color-flavor locked superfluid. At asymptotically high density, and thus weak 
coupling limit, all the coefficients of the phonon effective field theory can be computed from QCD. A studyof the 
shear viscosity and of the bulk viscosity coefficient C2 due to phonons has been already presented in Refs. [20l . [2l| . 
The contribution of kaons to C2 has been studied in Ref. [2^ . [23j . The contribution to thermal conductivity due to 
phonons was first estimated in Ref. [l^l , and recently computed in Ref. [2^ . 

This paper is organized as follows. In Section[TT]we review the hydrodynamical equations for a relativistic superfluid, 
when effects of dissipation are included. In Section [1111 we present the effective field theory for the superfluid phonon, 
at the leading and next-to-leading order, following the work of Ref. . In Section lTVl we derive the expressions for 

the phonon contribution to the bulk viscosity coefficients. Based on these formulas, and on the scattering amplitudes 
for the pertinent processes, we present a general discussion on the scaling behavior of these transport coefficients in 
Section |Vl Section IVII is devoted to the study of the CFL superfluid. We present our conclusions in Section IVIII 
We leave for the Appendix the proof that one needs the phonon dispersion law beyond linear order to obtain the 
first non-trivial contribution to the bulk viscosities. Throughout we use natural units, h = c = ks = 1, and metric 
conventions — , — , — ). 



II. HYDRODYNAMICS OF A RELATIVISTIC SUPERFLUID 

The hydrodynamical equations of a relativistic superfluid have been derived using different formulations 
Here we will use the one derived by Son [T3 |. The superfluid properties of a system arise from the spontaneous 
breaking of a continuous U{1) symmetry, with the appearance of a Goldstone mode. Since hydrodynamics is an 
effective field theory valid at long time and long length scales, the standard fluid variables couple to the Goldstone 
field. 

The hydrodynamical equations for the superfluid take the form of conservation laws for both the current, n^, and 
energy- momentum tensor, T'^", of the system 

9^71^ = , df^T^" = . (1) 

One further adds the Josephson equation, which describes the dynamical evolution of the Goldstone fleld, tp, or phase 
of the condensate 

u'^df.ip + M = , (2) 

where /i is the chemical potential of the system (26j . 

In this formulation of the superfluid hydrodynamics there is a clear interpretation of both the current and the 
stress-energy tensor as being due to the sum of the normal fluid part and the coherent motion of the condensate (the 
superfluid), as these are expressed as 

T^"" = (p + P) u^'u" - Pg^"" + V^d^'ipd'^ip , (3) 
^ n^u" - V'^d''(p , (4) 
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where m'^ is the hydrodynamical velocity, p stands for the energy density, P is the pressure, and V is a variable 
proportional to the quantum condensate. The energy density obeys the relation p = ST + n^p — P, where S is the 
entropy of the system. 

The dissipative terms associated to the above hydrodynamical equations have been constructed in [l^, showing 
that in the non-relativistic limit they correspond to those appearing in Landau's two-fluid model. To this end it is 
better to write the hydrodynamical equations in terms of the new variable 

u;^ = - {d^ip + pu^') . (5) 

Then, one can show that in the non-relativistic limit the spatial component of this four vector becomes the counterflow 
velocity of Landau's equations, w = to(vs — Vn), with the identification Vs = [Tsj . 

In order to construct the dissipative terms in the fluid equations one has to define a comoving frame. One possible 
choice is that in the frame where ~ (1, 0, 0, 0) the particle current and energy arc fixed as 

n^' ^ {n, V^Vip) , T"" = e . (6) 

Now one imposes some restrictions to the possible dissipative terms that might be added to the current and energy- 
momentum tensor, and T^'', respectively. In particular, one can choose constraints similar to those of the Eckart 
frame for the normal fluid hydrodynamics, = 0, and u^u^T^'^ = . 

Dissipative terms are taken into account as follows. First, one modifies the Josephson equation 

w^a^^ ^-p-x, (7) 

with 

X = -C^^^,ly^w^') - dd^u'^ . (8) 
One also adds dissipative terms to the energy-momentum tensor, which take the form 

T^" = K {A'^-'u'' + A'^-'u") {d^T + Tu^dsu^) + rjA''^ A""^ (^dsu^ + d^ug + ^g^sdc^u"^ (9) 

where In principle, more terms respectful with the symmetries of the problem arc possible, 

but those are neglected, assuming that they are small. One simply retains six transport coefficients, namely the 
shear viscosity 77, the thermal conductivity k, and CiiC2jC3iC4 which are the four bulk viscosity coefficients. In the 
non-relativistic limit the expressions of the bulk viscosity coefficients agree with those introduced by Khalatnikov 01 , 
up to some mass factors 

According to the Onsager symmetry principle [H, the transport coefficients satisfy the relation (i = C4j while 

the requirement of positive entropy production imposes that ^,77, C2,C3 ^-re positive and that C^f < C2C3- 

The friction forces due to bulk viscosities can be understood as drops, with respect to their equilibrium values, in 
the main driving forces acting on the normal and supcrfiuid components. These forces are given by the gradients of 
the pressure P and the chemical potential p. One can write in the comoving frame 

P = Pcq-Cidiv(yV)-C2divu, (10) 
M = Moq - C3div(y^w) - Cidivu, (11) 

where Peq and Peq are the equilibrium pressure and chemical potential. Note that the drops in the driving forces 
proportional to Ci and (^3 lead to entropy production only when divw ^ 0. Indeed in case divw = 0, the entropy 
production rate is given by 0, ITEj l 

i? = C2(divu)2. (12) 

Therefore, when divw = 0, in the Josephson equation ([7]) one has that x — -^CidjiU^, but such a term does not lead 
to dissipation. 



III. EFFECTIVE FIELD THEORY FOR THE PHONON OF A RELATIVISTIC SUPERFLUID 

All the phonon properties that are needed for the evaluation of the bulk viscosities can be extracted from the low 
energy effective fleld theory associated to the rclativistic supcrfluid. This has been constructed in Refs. [iH[3j and 
we review it here. 
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The effective field theory for the Goldstone field is constructed as a power expansion over derivatives and over 
fields, ~ allowing only the terms that respect the underlying symmetries of the system one is considering. 

The coefficients appearing in the low energy Lagrangian can be in principle computed from the microscopic theory, 
through a standard matching procedure. 

Let us review the form of the effective field theory. Consider a relativistic system that experiences the spontaneous 
breaking of a U{1) symmetry. If 6 is the phase of the field that gets an expectation value in its ground state 
^ = |^|e^*®, then 8 is related to the Goldstone field ip by the relation Q = iJ-ot — tp, where fiQ is the zero temperature 
chemical potential [l^, and t is the time. 

At the lowest order in 0, the Lagrangian must be a function of 

X - Ig'^'d^ed.e , (13) 

where g'^'^ is the metric tensor. Since we are not interested in curved space-time we can rewrite 

X = ^{ljl-2^iodop + d,ipd''ip) . (14) 
The Lagrangian for the Goldstone field at lowest order (LO) is expressed as 

£lo = P[V2X] , (15) 

where P is the zero-temperature pressure of the system one is considering. The reason why the Lagrangian takes this 
form is that the effective action of the theory at its minimum for constant classical field configurations has to be equal 
to the pressure. 

One can expand the pressure function around /ig finding in this way the Lagrangian 

^Lo = lido^? - y (5.^)2 - gsdopid^pd-'p) + 9i{d,pd''(pf + • • • , (16) 

where we have rcscaled the Goldstone field so as to have a kinetic term properly normalized. We have also neglected 
a piece proportional to d^Lp, only relevant for vortex configurations. The constants 0^,33,(74 can be written in terms 
of different ratios of derivatives of the pressure. In particular, is seen to agree with the zero temperature speed of 
sound of the system 

"'^7, {w) {m) ' ^^^^ 

while the coupling constants are expressed as 

_j_(d^ I dp_\ 1 _j_f^ 1 dp\ fd^py^ 

At next-to-leading order (NLO) there are four different structures, but two of them vanish in flat space-time. For 
simplicity, we will not write them down. At NLO one has 



^NLO = ^^d.Xd'^X + h{X){d,d''Qf , (19) 
Mo 

and we have re-defined the functions that appear in Ref. [Tgj so as to make /i and /2 dimensionless functions. 

We will mainly be concerned with the effect that the NLO terms have in modifying the phonon dispersion law. To 
this end, it is enough to expand the Lagrangian up to terms which are quadratic in ip. We assume that the functions 
/i^2 that appear in Eq. ([19]) are analytical, and that they can be expanded in a Taylor series around Hq. Thus, after 
the necessary re-scaling to have a normalized kinetic term, one finds 

-Clo + >Cnlo = lidopf - jid.ipf + mi {d^d^^d^doip) + m2{d,d''^f + ■■■ (20) 



where 



"^l = Q2p , W2 = gap , (21) 
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which leads to the following modification of the dispersion equation for the phonon 

- c^q^ + 2TOia;2(w2 _ ^2^) _^ 2m2(a;2 - = . (22) 

Assuming that uj,q<^ . / then one finds the dispersion law 

V ^0 



e, = Csq + Bq^ + ■■■ (23) 



where 



B = (^c,(l - cl)mi - 1(1 - c2)2m2^ . (24) 

The NLO Lagrangian also introduces corrections to the phonon self-couplings, that we will not explicitly write 
down here as they will only be necessary if we were interested in computing corrections to the leading order behavior 
of the transport coefficients. 

It is also interesting to discuss the conformal invariance of the phonon effective field theory. If the macroscopic 
system is conformal invariant, this symmetry should be inherited by the low energy effective field theory. At LO, 
this information is encoded in the explicit form of the pressure of the system, and correspondingly, in the speed of 
sound that should take the value of l/-\/3- At NLO, conformal invariance puts some constraints on the form of the 
functions /i and /2, which can only depend on powers of the ratio X / If conformal invariance is not a symmetry 
of the system, the functions P, /i and /2 will depend on the scale responsible for the breaking of that symmetry. 

It is important to stress the general character of the present discussion, valid for any relativistic superfluid, indepen- 
dent of the fact that the system is strongly coupled. Indeed in the above derivation one only employs the symmetries 
of the system. As an application one may think to the neutron superfluid which is believed to exist in the interior of 
neutron stars, or to the CFL superfluid that we will treat more carefully in the last part of the present manuscript. 

IV. PHONON CONTRIBUTION TO THE FREQUENCY DEPENDENT BULK VISCOSITIES 

At very low temperatures phonons give the leading thermal contribution to all the thermodynamic properties of 
the superfluid. In the hydrodynamic regime, phonons also give the leading contribution to the transport coefficients 
entering into the two-fluid equations. Khalatnikov developed the kinetic theory associated to these degrees of freedom 
for a non-relativistic superfluid ■ One can also construct the transport theory associated to these Goldstone mode 
excitations of a relativistic superfluid (27l . [29j . for example using a gravity analogue model, showing that in the 
non-relativistic limit one recovers the theory of Khalatnikov [23] . 

For the computation of the bulk viscosity coefficients we will use the method suggested by Khalantnikov for su- 



perfluid "'He, and recently used for the cold Fermi gas in the unitarity limit 17[. In this approach one studies the 
dynamical evolution of the phonon number density A/"nVi, and relates the variation of this quantity with the transport 
coefficients of the system. It has been shown in Ref. [I^l, that this approach is equivalent to solving the phonon 
transport equation in the relaxation time approximation. Note that with the RTA one obtains the correct parametric 
behavior of the transport coefficients, but one cannot fix with precision the numerical factor in front of these quan- 
tities. In order to obtain such a factor, it is necessary to solve numerically the corresponding Boltzmann equation, a 
rather complicated task. For the astrophysical applications we have in mind, getting the correct parametric behavior 
of all the transport coefficients is good enough, as this fixes the time scales of different physical processes. 

Consider a perturbation to the system that changes the number of phonons per unit volume A/'ph- CoUisional 
processes will tend to restore the equilibrium phonon number density. The dynamical evolution of A/'ph can be extracted 



from the transport equation obeyed by the phonon distribution function |27| . If we further use the relaxation time 



approximation to deal with the collision term, and assume that the speed of sound is constant, then we obtain that 

a,(K''A/-ph) = -^ , (25) 

Trol 

where 5M-p\i measures the departure of the phonon density from its equilibrium value. The relaxation time is given by 

— (26) 



where Fph is the rate of change of the number of phonons, to be specified later on. 
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From now on we will work in the comoving frame = (1,0,0,0) to simplify the treatment. In equilibrium, the 
number density of phonons is a function of the total current density, n = uq + V^hq, and of the entropy S. Since 
we are interested in astrophysical applications we shall assume that the perturbation is periodic, so that there is a 
dependence on time of the sort ~ exp {iiJct), with typical frequencies Wc, as measured in the comoving frame. For 
astrophysical applications Wc is of the order of the frequency of rotation of the star. 

Since the phonon number density depends on the current density and on the entropy we have that 



SK,^^Sn+^SS. (27) 
on ob 



The hydrodynamic equations take the form 



uf^d^n = -nd^^u^ - d^iV^Lj'') , (28) 
u^'d^S = -Sd^u^ , (29) 



and therefore in the comoving frame one can rewrite the phonon number variation as 



<5AAph - AAph - - divu - ^div(FV) . (30) 

At this point one studies how the variation of the phonon density alters the equilibrium values of the pressure and 
chemical potential Q . In this way one can identify the different bulk viscosity coefficients and if we define 



we obtain that 



Ci = , ^hh , (32) 

C2 = . T^/l , (33) 

1 - lUJcTrcl Aph 

C3 = (34) 

1 - iWcTrcl Aph 

and therefore (,1 = C2C31 meaning that one of the relation for positive entropy production is saturated. 

For Wc = and in the non-relativistic limit we recover the expressions of the bulk viscosity coefficients that 
were obtained by Khalatnikov For non- vanishing frequencies we obtain bulk viscosities that have both real and 
imaginary parts, as it occurs when there are periodic perturbations in time in the system, see [l[. For the computation 
of the dissipated energy in the system, only the real part matters (see for example [1^). One can see that the various 
coefficients attain their maximum value when a;^ = 

At this point it is easy to show that for phonons with a linear dispersion law all the bulk viscosity coefficients 
vanish, independent of whether the system is conformal invariant or not, see the Appendix for a short proof. This 
is in agreement with the results obtained for ^He [soj , and for cold Fermi atoms close to the unitary limit [l7| . Then 
we consider the first correction to the phonon dispersion law given by the cubic term, in the form given by Eq. \2?>\ . 
After defining the dimensionless parameter 

X = , (35) 

we compute the various quantities evaluated with the equilibrium phonon distribution function to the leading order 
in X. In this way we obtain that the number of phonons per unit volume is given by 

Kh = ^ (r(3)C(3) - X r(6)C(5) + Oix')) , (36) 



while the entropy turns out to be 



5. 



T 



3 



ph 



(r(5)C(4)-3xr(7)C(6)+0(x2)), (37) 
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where r{z) and ({z) stand for the Gamma and Riemann zeta functions, respectively. 

For the evaluation of the functional derivatives that appear in Eq. ([51]) we use the same strateg y fo llowed in Ref. ■ 
We use as independent variables the temperature T and the zero temperature chemical potential [2a| . Using a Maxwell 
relation, and the phonon contribution to the chemical potential, we arrive at the expressions (see Ref. [l7| for more 
details) 

and 

For a conformal invariant non-rclativistic superfluid one expects that Ci and (2 vanish, while there is no constraint 



on |3l| . For the ultracold Fermi superfluid in the unitarity limit this is what one indeed finds |17| . One can expect 
that these results hold for relativistic superfluids, as well. In the following we shall explicitly prove that this is what 
happens. 

Let us study under which conditions in a conformally invariant system the phonon contribution to the first and 
second bulk viscosity coefficients vanish. From the expressions in Eqs. ((33)) and (|34p this is equivalent to asking that 
I2 = 0. In a conformally invariant system the speed of sound is fixed, = 1/3. Further, since /xq is the only energy 
scale it follows that by dimensional analysis B oc Cs//Xq and n ~ /Xq. Upon considering these scaling behaviors in 
Eq. ([M)) . and in the case that B ~ uCs/^Iq one finds that I2 vanishes. 

In Scction lVTl wc will discuss the CFL superfluid, where we find that for masslcss quarks, and neglecting the running 
of the gauge coupling constant, the parameter B has the dependence expected in a conformal invariant system, and 
thus the only non- vanishing bulk viscosity turns out to be (.3. 



V. ZERO-FREQUENCY TEMPERATURE DEPENDENCE OF RELAXATION TIME AND BULK 
VISCOSITIES OF A COLD RELATIVISTIC SUPERFLUID 

In the previous Section we have derived the expressions of the three independent bulk viscosity coefficients. In order 
to determine their actual temperature dependence we now estimate the temperature dependence of the relaxation 
time for the pertinent process. We shall consider the situation where Wc = 0; the non-static case can be easily deduced 
as well. For the time being the only relevant assumption we will make is that the contribution of other quasiparticles 
to transport phenomena is thermally suppressed. Therefore, we shall restrict to the situation where T is much less 
that any energy gap of the system. 

The computation of the phonon contribution to the bulk viscosities presents different subtleties. As we saw in 
Section lllll the phonon self- interactions are described by a low energy effective field theory, whose form is the same 
for all the systems that share the same global symmetries. However, the coefficients that appear in the low energy 
Lagrangian depend on the system considered, as their values have to be matched with the corresponding microscopic 
theory. As we saw, the bulk viscosities turn out to depend on the coefficient B that appears in the phonon dispersion 
law, see Eqs. ([55)) and ([55)) . The transport coefficients also depend on the relaxation time of the coUisional process 
responsible of the corresponding dissipative phenomena. For bulk viscosity those are collisions that change the number 
of phonons. It turns out that the sign of B decides whether some collisional processes are kinematically allowed or 
not, and ultimately, this affects how the relaxation time associated to bulk viscosity scales with T. 

For a dispersion law where B is positive, the so-called Beliaev process [13 which describes the decay of one phonon 
into two, if — > ifip, is kinematically allowed. The decay rate associated to this process is given by 

Fph- J dPdA-dQ|Al(P,g,/C)|Vq(po)(l + ricq(go))(l + rieq(fco))(27r)4<5(4)(p_ ^ (49) 

where = {po,p) is the four momentum and we have defined the integration measure as 

JdP^I ^2Q(j>,)5{pl-el), (41) 

and rieq is the Bose-Einstein equilibrium distribution function. The squared of the scattering amplitude at LO is 
expressed as 



\M{P, Q, K)\^ = gl {po Q-K + q„P- K + koQ- Pf , 



(42) 
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written in terms of the momenta of the phonons and of the three phonons self-couphng constant defined in Eq. (|18|) . 
Even without an exphcit evaluation of the above integral, one can infer that Fph ^ g^T^/c^, and that the relaxation 
time scales with the temperature as Tioi ^ cl/g^T^. 

As we have already shown in the previous Section, for a conformally invariant system, Ci = C2 = Oj whereas C3 7^ 0. 
Upon substituting the scaling behavior of Trci, of /i and I2 given in Eq. (p8| and ([39)) . of J\fph given in Eq. ((36)) . in 
Eqs. ([34| one finds that C3 ^ T^. When conformal symmetry is broken, it follows that the three independent bulk 
viscosity coefficients do not vanish and have the same temperature dependence, i.e. Ci ^ T^. 

Including NLO corrections both in the dispersion law gp, and in the self-coupling 33, corrects this leading order 
behavior with terms of order B^T^/c^, which at low temperatures are negligible. 




FIG. 1: One of the Feynman digrams for the process ipLp — )■ ipip^p. The scattering matrix of this process is reported in Eq. (|44|) . 



For a phonon dispersion law that curves downward, B < 0, the Beliaev process is not kinematically allowed, and 
one has then to consider the five phonons coUisional processes. The decay rate for these processes is very different 
with respect to the </? — ipip process. Further, as we discuss below, collisions with large or small scattering angle also 
have a rather different decay rate. 

The decay rate of the process ip(p — >■ (pip(p is given by 

rph= / n dP^M5^\^n,^{pl)n,^(jjl){l + n,^ipl)){l + n,^(j^^^^^ 

•' i=l,..,5 

(43) 

where is the scattering matrix for the process considered. 

After a careful analysis, one realizes that this rate has a very different behavior if one considers large or small angle 
collisions. Let us consider the Feynman diagram reported in Fig. [TJ The square of the scattering matrix behaves as 

lA/f (p p p p p _ 2 2 PiiPi,-P3,-Q)FiiP2,Q,-P4,'-P5) f... 

\J^5:p[Pl,P2,P3,P4:,P5)\ -.9354 — ^72 ' I44j 



•,2 _ g2 



where Q'^ = (^OiQ) ~ (Pio ^PsoiPi ~ Ps) is the momentum transfer, and and F4 are polynomial of the phonon 
momenta that we shall not specify here, but that can be easily inferred from the Lagrangian (fT6| . If one considers 
a collision with large angle scattering, naive power counting suggests that Fph ~ glgfT^^ /cl'^ , and thus t,.oi ^ 
c^/{g^glT^^), which turns out to be a rather large value for low temperatures. 

Processes associated to small angle scattering behave very differently. Let us start by considering phonons with a 
linear dispersion law. Then, the scattering matrix associated to a small angle collision behaves as 

lA/f (p p p p p \\^ - 2 2 Pi{Pi^-Pi^~Q)FliP2.Q,~Pi.~P5) 

\J^bv\Pl, P2, Pi, Pb)\ - g?,gA /I ^4^2^2/1 \2 ' v45J 



4cfp2p2(i„cos6'i3j 

where Oiz is the angle between the vectors pi and P3. In the limit where ^13 — > this process presents a collinear 
singularity, which needs to be regularized. This might be done by considering a phonon dispersion law that includes 
a cubic term, as then in the denominator of Eq. (j45p one has a term proportional to i?^, which does not vanish in the 
013 limit. 

Let us stress here that while we have based our discussion on the Feynman diagram in Fig.[Tl there are other channels 
for five phonons collisions that also present similar collinear singularities. The naive power counting associated to small 
angle collisions then suggests that Fph ^ g^g^T^^ / (c^B^), and thus a relaxation time scaling as TioI ^ c^B^ / (g^g^T^), 



9 



which is a much faster process than that corresponding to large angle collisions. Because in the computation of the 
transport coefficients one has to consider first the collisions with the shorter relaxation time, one then concludes that 
in this case, and when conformal symmetry is broken, all the bulk viscosity coefficients scale with the temperature as 
Ci ^ y. Corrections to this leading behavior would be very suppressed, as discussed earlier. 

Summarizing, bulk viscosities in the cold regime of relativistic superfluid strongly depend on the form of the phonon 
dispersion law. If conformal symmetry is unbroken, one should have Ci = C2 = 0; while the three bulk viscosity 
coefficients are non-zero when conformal symmetry is broken. When B > the process ip — > ipip is kinematically 
allowed and gives the dominant contribution to bulk viscosity. In this case one may expect a rather small value for the 
bulk viscosities, with scaling Q ^ T'^. When B < one has to consider the processes (p(p -> ipipip which is dominated 
by small angle scatterings. These processes give scalings d ^ 1/T. Let us finally remark that while our discussion 
has been done for relativistic superfluids, it is also valid for non-relativistic superfluids. 



VI. BULK VISCOSITIES FOR THE CFL PHASE DUE TO SUPERFLUID PHONONS 

At asymptotically high densities quark matter is believed to be in the color-flavor locked phase In this phase 
up, down and strange quarks of all three colors form zero-momentum spinlcss Cooper pairs. The CFL pairing 
pattern spontaneously breaks both the flavor chiral and the baryonic number U{1)b symmetries, leaving a discrete 
Z2 symmetry. There arc eight (pseudo) Goldstone modes associated to the chiral group, and one exact massless 
Goldstone boson associated to the breaking of the U{1)b, which is responsible for the superfluid properties of the 
system. All gluons are massive, and all quarks are gapped. All the properties of CFL quark matter can be computed 
from QCD in the high density limit 

In this Section we derive the bulk viscosity coefficients of CFL quark matter using the strategy discussed in the 
first part of the manuscript. 



A. Conformal limit 



As shown in Section IIIIl the effective low energy Lagrangian for the Goldstone boson can be obtained from the 
pressure and for the CFL quark matter at asymptotic densities we have that 

P[Mo] = , (46) 

where hq is the quark chemical potential. At very high /ig, where the coupling constant is small .g(/io) <C I, the effects 
of interactions and the effects of Cooper pairing are subleading and have been neglected in Eq. (|46| . Also, at very 
high densities one might neglect the quark masses, as niq fiQ. 

The values of all the coefficients which appear in the phonon Lagrangian at LO, Eq. (|16p . can now be extracted as 
explained in Section [llll and are given by 

It is possible to compute the phonon dispersion law at NLO for massless quarks by matching with QCD at finite 
density [s^ [s^l ■ In this way one obtains that 

where A is the CFL gap. In the asymptotic high density limit, the gap can be computed from QCD [s^ 

A ~ fxog-'' exp , (49) 

where g is the QCD gauge coupling constant, and 60 = 5127r'*(|)^/^ cxp [svllssj. 

For sufficiently large chemical potentials, one can neglect the running of the coupling constant [s^. Then one can 
consider that the gap A oc ^o- that is, the gap only depends linearly on the chemical potential and therefore docs not 
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break conformal invariance. Moreover in the regime we consider the total particle density is n w consistent with 
the approximation to the CFL pressure, and being in the regime T <C /io- Since n oc /Iq and A cx /io, it follows that 

/i = -^T^(-^C(3)-7C(5))^, 1,^0, (50) 

where, according with Eq. (j48p . the coefficient of the cubic term in momentum is given by 

Since _B < 0. we have to consider that the processes giving the leading contributions to the bulk viscosity are small 
angle five phonons collisions. The estimate of the relaxation time for these processes is given in Sec. IVland turns out 
to be 

„3 ,,12 

We consider first the case Wc = 0. From Eqs. ((50)) we conclude that Ci = C2 = 0, while the third bulk viscosity does 
not vanish and depends parametrically on the physical scales of the problem as 

We remark that these arc only approximated results that arises in the 5 <C 1 limit, after neglecting the running of the 
QCD gauge coupling constant and the effect of the strange quark mass. These are good approximations in the regime 
of high density [39| . For non- vanishing values of oJc one can accordingly find the value of the only non- vanishing 
transport coefficient, given the expressions of the relaxation time and zero- frequency viscosity, see Eq. p4p . 



B. Including scale breaking effects due to the strange quark mass 

The bulk viscosities are sensitive to scale breaking effects. The quantum scale anomaly breaks the conformal 
symmetry introducing, through dimensional transmutation, the quantum scale Aqcd- One could then compute g- 
corrections to the pressure of quark matter, which would then modify the different coefficients of the phonon effective 
field theory, introducing terms proportional to the QCD beta function. In the very high /xq limit, when g{fio) ^ 1, we 
expect this to be a rather negligible effect, as mentioned in the previous Section. Also consider that for large values 
of the chemical potential the gap parameter. A, does not break scale invariance because in this case one has A cx /xq. 

The inclusion of quark mass effects leads to the breaking of scale invariance. In the CFL phase, the three light 
quarks participate in the pairing process, and therefore the largest effect comes from the strange quark mass m^. 
However, one cannot consider arbitrarily large values of the strange quark mass, because otherwise the CFL phase 
becomes chromo-magnetically unstable [43[. Therefore, we shall take ml < 2A^0i which is the threshold value under 
which the CFL phase is stable. Since ^ /ig we will consider the leading order corrections in m^Z/XQ. After imposing 
the constraints of electrical neutrality and beta equilibrium of quark matter, the first correction of the order m^/^o 
to the pressure reads Q 

^[a^o] = ^ (/"o - Mo^s) • (54) 

Using Eqs. ([T7| and (fTS]) we find that to the order ml/ ji^ the coefficients of the LO phonon Lagrangian are corrected 
due to this scale breaking effect as 

f TTl^ \ f Tfl^ \ f m^ \ 

c^=cL(^l-5^j , .3=,93,o(^l + ^j , 9. = 9..o[l + ^J, (55) 

where the Cs^, 173,0 and (74.0 refer to the values obtained when m^ — 0. Thus, the leading order effect of the strange 
quark mass in the phonon Lagrangian is to modify the velocity of the phonon, that is the speed of sound, and a finite 
renormalization of the cubic and quartic phonon self-couplings. 

The phonon dispersion law at cubic order including quark mass corrections has not yet been computed in the 
literature, as it is subtle. Because a mass term in the QCD Lagrangian couples the particle and antiparticlc degrees of 
freedom, that computation requires the knowledge of the antiquark propagator in the CFL phase, exactly as it occurs 
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in the computation of the CFL meson masses [40|-|42[. In any case, the correction of the coefficient of the cubic term 
in momentum in the phonon dispersion law will be of the form 

B = Bo + , (56) 

with bs a dimensionless constant. Since we are only interested in determining the scaling laws and not the numerical 
factors, we leave the computation of this constant for a future publication. 

Upon substituting the expression of Eq. ([55]) and (|56)) in Eqs. (|38)) and (|38|) . one finds 

20 /5 



f^T^ (tt^CCS) - 7C(5)) 2Boc.,o (| " ) ^ ■ (58) 



7cV2 V ' 'J --vs 7 Mg 

The scale breaking terms affect the relaxation times as well, because the phonon self-couplings are modified for 
TOg 7^ 0, as shown in Eqs. (|55p. A naive power counting analysis suggests that 



T-..ci^rreLo(l + (25,,-i)^) . (59) 



We discuss now the values of the bulk viscosities for the Wc = case, the non-static case is then easily deduced as 
well. As we are only interested in the scaling laws associated to the bulk viscosity coefficients, we conclude that to 
order ml/fif, 



1 Mo rn^s . 1 Mo^ rnj 



Ci-^T^^^ C^-^T^"^- (60) 



Further, taking into account that A oc /io, one obtains the dependence of the different coefficients with the basic 

2 4 

physical scales of the problem, T, /io,™^- In this way one finds [Ci] = [C2] = ^ ^nd [C3] — Note that 

the scaling behavior of C2 was previously found in Ref. j2l| . There the phonon Boltzmann equation was solved 
numerically, using a variational method. In that article it was assumed a phonon spectrum linear in momentum, and 
the interactions were extracted from the LO Lagrangian only, neglecting all phonon physics at NLO. In the linearization 
of the transport equation, the variations of the speed of the phonon with the total density were neglected, and that 
is why a non-zero result was obtained, even if the whole computation was carried out with the LO phonon physics. 



VII. CONCLUSIONS 



Superfiuidity arises at low temperatures after the spontaneous breaking of a global U{1) symmetry. It is possible 
to determine the expression of the low energy Lagrangian of the Goldstone mode by using symmetry arguments. 
However, the coefficients of the effective field theory are not fixed by the symmetries of the system; they have to be 
matched with the microscopic field theory, and thus depend on the physics at the very short scales. 

In this paper we have used the effective field theory describing the low energy properties of a superfluid system, 
together with kinetic theory in the relaxation time approximation, to obtain the contribution of the superfiuid phonon 
to the bulk viscosity coefficients. We have first presented a general discussion on the low energy Lagrangian of the 
system. Then, we have determined the scaling behavior of the bulk viscosities with T. In order to do this we have 
estimated the relaxation time of the processes responsible for bulk viscosity. The allowed processes depend on the 
sign of the correction to the linear dispersion law of the phonon. If the correction to the linear dispersion law curves 
upward, one phonon can decay into two, but this process is kinematically forbidden in the opposite case. If the 
correction to the linear dispersion laws curves downward, then the most important process is iptp — > (fifip which is 
dominated by small angle scatterings. We have presented a general discussion on how the relaxation times behave in 
both cases in Section IVl 

In the derivation of the transport coefficients we have neglected the contribution of other quasiparticles. This 
approximation is justified in the situation where T is much less that any other energy gap of the system. If not, one 
should take into account the contribution of other quasiparticles as well. 

We have used the general expressions of the transport coefficients for a relativistic superfluid to determine the 
phonon contribution to the bulk viscosity coefficients of CFL quark matter. In order to do this we have determined 
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the scaling of the relevant relaxation times. Let us stress that with the relaxation time approximation we can only 
obtain the correct parametric behavior of the transport coefficients, but this should be enough to get their correct 
scales to be used, for example, in the study of the r-mode oscillations of a quark star. Previous analysis for the 
r- modes of a CFL quark star have been carried out in Refs. j43 - l46| . but there the transport coefficients that we have 
studied as well as the mutual friction that operates in rotating superfluids [47j were ignored. Let us also point out 
that the knowledge of the bulk viscosities is relevant for the evaluation of the cavity conditions in quark stars [i^ . 

It would be interesting to study the phonon contribution to the bulk viscosities of a neutron superfluid applying the 
general expressions derived in the present paper, as that contribution has been ignored in astrophysical applications. 
It might also be interesting to study the phonon contribution to transport phenomena with holographic techniques, 
as this should be independent of whether the underlying microscopic theory is strongly or weakly coupled. 
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Appendix A: Phonons with a linear dispersion law 



For phonons with a linear dispersion law one can write 

_ r(3)c(3) 

meaning that A/'ph = J^phiS) and the phonon density is independent of n at this order. From Eqs. p2l33l34p it follows 
that for phonons with a linear dispersion law Ci = C2 = Cs = 0- Note that this result is true for any form of the speed 
of sound, Cs = Cs(/io). The fact that for phonons with a linear dispersion law A/'ph is independent of n can also be 
explicitly proven. By the chain rule we have that 

dn dT dn dfiQ dn ' 

that we can rewrite as 

1 dAf^_^fdT\ ^fdcs^ ^^3^ 



(A4) 



3Afph dn T \dn J g Cg \dn 
From the value of the entropy with a linear phonon dispersion law we have that 



/ r(5)c(4) y/^ T 



therefore 

Cs \dn J g T \dh 

and it follows that = 0. 

on 



(A5) 



References 



[1] L. Landau and Lifschitz, "Fluid Mechanics" vol. 6, Prentince Hall, New Jersey. 

[2] L M. Khalatnikov, "Introduction to the Theory of Superfluidity", Benjamin, New York, 1965. 



13 



arXiv: astro-ph/0 1 01136" 
arXiv;0803.3295 [hep-th]] 



[3] For reviews, see K. Rajagopal and F. Wilczok, arXiv:hep-ph/0011333 ; M. G. Alford, Ann. Rev. Nucl. Part. Sci. 51, 131 
(2001); G. NarduUi, Ri v. Nuovo Ci m. 25N3, 1 (2002); M. G. Alford, A. Schmitt, K. Rajagopal and T. Schafer, Rev. Mod. 
Phys. 80, 1455 (2008) [arXiv:0709. 4635 [hep-ph]]. 

M. G. Alford, K. Rajagopal and F. Wilczek, Nucl. Phys. B 537, 443 (1999) [arXiv:liep-pli /9804403 . 
P. W. Anderson and N. Itoh, Nature 256, 25 (1975). 

N. Andersson and K. D. Kokkotas, Int. J. Mod. Phys. D 10, 381 (2001); L. Lindblom, 
S. A. HartnoU, C. P. Herzog and G. T. Hor owitz, Phys. R ev. Lett. 101, 031601 (2008) 
C. P. Herzog, P. K. Kovtun and D. T. Son, larXiv:0809.487 [hep-th]. 

V. V. Lebedev and I. M. Khalatnikov, Zh. Eksp. Teor. Fiz. 56, 1601 (1982); I. M. Khalatnikov and V. V. Lebedev, Phys. 
Lett. 91A, 70 (1982). 

B. Carter and I. M. Khalatnikov, Phys. Rev. D 45, 4536 (1992). 
Carter and D. Langiois, Phys. Rev. D 51, 5855 (1995). 

Carter and I. M. Khalatnikov, Annals of Physics 219, 243 ( 1992). 

Andersson and G. L. Comer, Living Rev. Rel. 10, 1 (2005) iarXiv:gr-qc/0605010 . 
J. Mod. Phys. A 16S1C, 1284 (2001) [arXiv:hep-ph/0011246|. 



T. Son, Int 

E. Gusakov, Phys. Rev. D 76, 083001 (2007) [arXiv:0704.1071l [astro-ph]] 
Onsager, Phys. Rev. 37, 405 (1931). 

A. Escobedo, M. MannareUi and C. Manuel, Phys. Rev. A 79, 063623 ('2009') [arXiv:0904.3023l [cond- mat. quant-gas], 

T. Son, arXiv:hep-ph/0204199, 

T. Son and M. Wingatc, Annals Phys. 321, 197 (2006) (arXiv:co nd-mat/0509786 . 



Manuel, A. Dobado and F. J. Llanes-Estrada, JHEP 0509, 076 (2005) ^arXiv:hep-ph/0406058 

Manuel and F. J. Llanes-Estrada, JCAP 0708, 001 (2007) [arXiv:0705 .39091 1 hep -ph] j . 

G. Alford, M. Braby S. Reddy and T. Schafer, Phys. Rev. C 75, 055209 (2007) 'arXivmucl-t h/0701067 
G. Alford, M. Braby and A. Schmitt, J. Phys. G 35, 025002 (2008) arXiv:0707.2389 [nucl-th]]~^ 



M. , „ 

I. A. Shovkovy and P. J. Ellis, Phys. Rev. C 66, 015802 (2002) ar Xiv:hep-ph/0204132| . 
M. Braby, J. Chao and T. Schaefer, arXiv:0909.4236 [hep-ph]. 

As explained in Ref. 0], Eq. (18-13), the chemical potential receives a contribution from the thermal excitations of the 
system. We shall refer to the chemical potential that does not include such excitations as the zero temperature chemical 
potential. 

M. MannareUi and C. Manuel, Phys. Rev. D 77, 103014 (2008) arXiv:0802.0321 [hep-ph]]. 
L. Lindblom and B. J. Owen, Phys. Rev. D 65, 063006 (2002) arXiv:astro-ph/0110558] . 
V. Popov, Gen. Rel. Grav. 38, 917 (2006) arXiv:gr-qc/0607023j: 
I. M. Khalatnikov and D. M. Chernikova, Sovie t Physics JETP 23, 274 ( 1966). 
D. T. Son, Phys. Rev. Lett. 98, 020604 (2007) [arXiv:cond-mat/051172T] . 
S. T. Beliaev, Sov. J. Phys. 7, 289 (1958); 34, 299 (1958). 
K. Zarembo, Phys. Rev. D 62, 054003 (2000) arXiv:hep-ph/0002123 



K. 
D. 

W 



The phonon dispersion law for the CFL superfluid phonon has also been computed numerically in Ref. [35|], but using as 
a microscopic theory a NJL model rather than QCD. In this case, the dispersion law shows the same behavior one has in 
QCD, that is B < 0. 

Fukushima and K. lida, Phys. Rev. D 71 , 074011 (2005 ) 'arXiv: hep-ph/0501276 
T. Son, Phys. Rev. D 59, 094019 (1999) 
E. Brown, J. T. Liu an 

Wang and D. H. Rischke, Phys. Rev. D 65, 054005 (2002) arXiv:nucl-th/0110016 
Schafer, arXiv:hep-ph/0402032, 



l arXiv:hep-ph /9812287 
c. Ren, Phys. Rev. D 61, 114012 (2000) 



arXiv:hep-ph /9908248 



T. Son and M. 
T. Son and M. 
Manuel and M. 
Casalbuoni, R. 



A. Stephanov, Phys. Rev. 
A. Stephanov, Phys. Rev. 
H. G. Tytgat, Phys. Lett 
Gatto, M. MannareUi, G 



arXiv:hep-ph/9910491 



D 61, 074012 (2000) 

D 62, 059902 (2000) 

B 479, 190 (2000) farXiv:hep-ph/0001095' 
NarduUi and M. Ruggieri, Phys. Lett. B 605, 



arXiv:hep-ph /0004095 



615, 297 (2005)] arXiv:hep-ph/0410401 . 
J. Madsen, Phys. Rev. L ett. 85, 10 (2000). 
B. A. Sa'd. larXiv:0806.3359 i [astro-ph]. 
P. Jaikumar, G. Rupak and A. W. Steiner, 
M. MannareUi, C. Manuel and B. A. Sa'd, 
J. Madsen. [axXiv:0909.3724i [astro-ph.SR]. 



362 (2005) [Erratum-ibid. B 



Phys. Rev. D 78, 123007 (2008) [arXiv:0806.1005] [nuc l-th]] . 
Phys. Rev. Lett. 101, 241101 (2008) [arXiv:0807.3264l [hep-ph]]. 



